
clear all
close all
%--------------------------------
  lw = 2.5; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 16;
  msize = 10;
%----------------------------------------
disp('From BUGS: alpha')

load alpha.txt
load beta.txt
load tau.txt
load sigma.txt
alph = alpha(:,2); bet = beta(:,2); ta = tau(:,2); sig = sigma(:,2);
figure(1)
histo(alph(alph < 6 & alph > 0),100,0,1)
axis tight
print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\alpha_bugs.eps'

figure(2)
histo(bet(bet < 3 & bet > -1),100,0,1)
axis tight
print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\beta_bugs.eps'
figure(3)
histo(ta(ta < 10),100,0,1)
axis tight
print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\tau_bugs.eps'
figure(4)
histo(sig(sig<5),100,0,1)
axis tight
print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\sigma_bugs.eps'
